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Abstract 

We investigate dual mechanisms for interacting particle systems. Generalizing an approach of 
Alkemper and Hutzenthaler in the case of coalescing duals, we show that a simple linear transforma- 
tion leads to a moment duality of suitably rescaled processes. More precisely, we show how dualities 
of interacting particle systems of the form H(A, B) = q^ AnB \ A, B C {0,1} N ,q 6 [—1,1), are rescaled 
to yield moment dualities of rescaled processes. We discuss in particular the case q = — 1, which 
explains why certain population models with balancing selection have an annihilating dual process. 
We also consider different values of q, and answer a question by Alkemper and Hutzenthaler. 

Keywords: Markov processes, duality, interacting particle systems, graphical representation, annihi- 
lation, selection. 

MSC Subject classification: 60K35. 

1 Introduction and main result 

Dualities have proved to be a powerful tool in the analysis of interacting particle systems and popu- 
lation models. For interacting particle systems, one generally considers two kind of duals: coalescing 
and annihilating duals, |Lig0 5 > Gri79l [SL95j . In connection with population models, rescaled interacting 
particle systems and their limits are of considerable interest, and it is natural to ask in which sense 
rescaling preserves dualities. Alkemper and Hutzenthaler |AH07j consider the case of coalescing dual 
mechanisms, and derive a 'prototype' moment duality under rescaling. Swart }Sw06j uses a similar idea 
to obtain dualities of stepping stone models. In this paper, we consider a general form of a duality for 
interacting particle systems, |SL95j . This includes coalescing as well as annihilating dual mechanisms. 
We prove a 'prototype' moment duality of linearly transformed rescaled processes, in a similar fashion as 
for the coalescing case, and we discuss the situation for annihilating duals in some more details. As an 
application, we consider branching annihilating processes and their duals. Our approach explains why 
population models with balancing selection generally have an annihilating dual process, as was found, for 
example, in jBEM07] . Since this moment duality is obtained directly from the graphical representation, 
our approach presents a pathwise construction. We finally discuss the connections to the Lloyd-Sudbury 
approach, [SL95 , Sw06 , which answers the question posed in |AH07j . 



For a Markov process (X t )t>o we write for the law of the process started in x, and E^, for 
the corresponding expectation. Two Markov processes {X t )t>o, (Y t )t>o with state spaces E and F, 
respectively, are called dual with respect to the duality function H : E x F — > R if for all t > 0, x £ E, 
y £ F the equality 

E x [H(X t ,y)]=E y [H(x,Y t )} (1) 
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holds. This means that the long-term behaviour of one process is - to some extent - determined by the 
long-term behaviour of the other process. The usefulness of a duality depends on the duality function 
H. If for example (X t ) takes values in R, and (Y t ) in N, we call a duality with respect to the function 

H{x, y) = x v 

a moment duality, since it determines all the moments of X t . For practical purposes, it is often useful 
to have a pathwise construction of the dual processes, which means that they can be constructed on 
the same probability space in some explicit way, for example as functions of one underlying driving 
process. In the case of interacting particle systems, this construction is usually provided by the graphical 
representation, IHar781 [Gri79j . We explain this below in the setup that we use for the present paper. 

Let N G N, and let En '■= {0, 1}^. We write x G Em as a vector x = {xi)\<i<N- A partial order on En 
is given by x < y Xi < yi VI < i < N. We write x Ay for the minimum of x and y with respect to this 
ordering. Let (X^) t >o and (Y t N ) t >o denote Markov processes defined on some probability space (f2, J- ', P) 
with values in E N , X t N = (X t N (i)\ =li ... iN . Let A? := {i : X t N {i) = 1} and 5f := {i : Y t N (i) = 1}; this 
defines processes taking values in the subsets of {1, ...,N}. We write \X X \ := l-A^I for the cardinality of 
the set A^ , that is for the number of l's. Sudbury and Lloyd SL95 argue that in this context, duality 
functions that are functions of A n B alone should be of the form 

H(A,B)=q\ AnB \, A,Bc{l,...,N}, 

for some q G R \ {1}. We take this as a motivation to say that two En— valued Markov processes 
(X t N ),(Y t N ) are q-dual if 

E x [q lx * AY « I] = E y [q\ x <> AY * I] Vx,y £E N , t> 0. (2) 

That is, the duality function is H{x,y) = q> xAy \. Special cases are q = 0, which is called coalescing 
duality, and q = — 1, which is called annihilating duality. In these cases the duality function becomes 
Ol^l = 1 {X A V =0}, and (-1)1^1 = 1 - 2 x l {[xAyl j s dd}' respectively. 

We now describe the graphical representation for such dualities. This is classical, for more details 
see |Har781 IGri79| . For each i G {1, N}, draw a vertical line of length T, which represents time up 
to a finite end point T. We consider ordered pairs (i,j) with i,j G {1, N}. For each such pair, run 
m G N independent Poisson processes with parameters (A^ ), k = 1, m. At the time of an arrival draw 
an arrow from the line corresponding to i to the line corresponding to j, marked with the index k of 
the process. Do this independently for each ordered pair, for each k = 1, ...,m. For each k, we define 
functions f k ,g k : {0, l} 2 -> {0,1} 2 . A Markov process (X^) with cadlag paths is then constructed by 

specifying an initial condition x = (xi)i—i n, and the following dynamics: X^ = x until the time of 

the first arrow in the graphical representation. If this arrow points from i to j and is labelled k, then 
the pair (xi,Xj) is changed to f k (xi,Xj), and the other coordinates remain unchanged. Go on until the 
next arrow, and proceed exactly in the same way. The dual process (Y t N ) is constructed using the same 
Poisson processes, but started at the final time T > 0, running time backwards, inverting the order of 
all arrows, and using the functions g k instead of f k . 

Following |AH07j . we call the functions f k ,g k basic mechanisms, and we generalize the definition 
of dual basic mechanisms given by Alkcmper and Hutzcnthalcr. For x = (xi,X2) G {0, l} 2 we use the 
notation x* := {x2,x\); the dagger accounts for the reversal of an arrow. 

Definition 1.1. Two basic mechanisms f, g : {0, l} 2 — > {0, l} 2 are called g-dual mechanisms if and only 
if 

q \xA(g(y^\ = q \f{x)Ay\ \/.T, y G {0, l} 2 . (3) 
It is clear that two processes constructed using q— dual mechanisms are q— dual processes. 
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Figure 1 



Lemma 1.2. Fir m e N, q 6 M \ {1} and T > 0. For euen/ fc = l,...,m, let f k ,g k be q-dual 
basic mechanisms. Let (X^), (Y t ) be Markov processes with state space En, constructed using the 
mechanisms fk, k = 1, . . . , m and gk, k = 1, . . . , m, respectively, and Poisson processes with the same 
symmetric parameters X k - = A^, k = 1, ....,m, i,j = 1, JV. TTiera i/iere exists a process (Y t N ) such that 



yN £ y n 



\X t AY \ _ 



g |XoAft| o.s. V0 < t < T. 



(4) 



Proof. Since we assume that the Poisson processes have the same rates, we can construct Y t from 
the graphical representation of (X^), using the same realization of the Poisson processes, reversing 



time and the directions of all the arrows. It is clear from the construction that then Y t N 



V N 
I t i 



and 



\X t /\Y \ 



-JXoAYtl hold ( see Figure 1). For some more details, in the case of coalescing mechanisms, 
compare the proof of Proposition 2.3 of jAH07j . □ 

Taking expectations, the following is then obvious. 

Corollary 1.3. In the situation of Lemma M.'A the processes (X^ 1 ) and (Y t ) are q—dual. 

Remark 1.4. We note that Lemma [1.21 tells us that (X^) and (Y t N ) are dual in a very strong sense, 
namely, the equation H(X^ , ig ) = H(Xq ,Y t ) holds almost surely instead of just in expectation. 
We call such processes strongly pathwise dual. We have just seen that a construction via graphical 
representation and q- dual basic mechanisms automatically leads to a strong pathwise duality. 

We are now ready to state and prove the main result of this article. 

Theorem 1.5. Let (X^), (Y t ) be Markov processes with state space Fjy that are q^ — dual for some 
ojv € [—1,1). Choose exchangeable initial conditions Xq ,1q £ En, fixing \X$\ = fcjy, \Yq\ = n^, and 
suppose that Xf^ and Y t N stay exchangeable for all t > 0. Assume that npj/N — > and E[|Y" t ^|/iV] — > 
as N ^ co, for some time scale tjy > 0. Then 



lim E 



N- 



1 + (q N - 1) 



\xl\ 

N 



\Y t N 



= lim E 

N->oo 



l + (q N - 1) 



N 



N 



provided that the limits exist. 

Theorem 11.51 applies, for example, to processes constructed from basic mechanisms with rates Ajj- 
that do not depend on i and j; all our later examples fall into this class. As illustrated by the Moran 
model, however, exchangeability of the arrows is not a necessary condition for the processes to stay 
exchangeable. 

Depending on the scaling, Theorem 11.51 may lead to a moment duality, if — >■ X t , and \Y t " \ — > Y t , 
as we then get E[(l + (q - l)X ) Yt ] = E[(l - (q - l)X t ) Ya ]. If X N and Y N have the same scaling, we 
may get a Laplace duality, see |AH07j for an example. 

Proof. The proof relies on the simple fact that the distribution of \X A Y\ given \X\ and \Y\ is approx- 
imately binomial with parameters \Y\ and provided that \Y\ is small with respect to N. Indeed, 
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\X A Y\ follows a hypergeometric distribution, since it is obtained by distributing the \Y\ l's of the Y- 
configuration onto the \X\ l's of the A-configuration, without hitting the same 1 twice. Approximating 
the hypergeometric distribution by a binomial distribution will give us the result. Let Z N ~ Bin (tin, ^y") 
with xn G {0, N} and n^/N — > 0. By Theorem 4 of |DF80j . we can bound the total variation distance 
between the hypergeometric and the binomial distribution as 



||Hyp(iV,a;Ar,n A r) - Bin(njv, -rf)\ 



N 



TV 



< 



4n N 
N ' 



Since we assumed qN G [—1,1), we obtain 



E 



In 



\X?J=x N ,\Y "\ 



TIN 



J2lNn\X t N AY \ = k | \X t N \=XN,\Y N \=n N ) 
=o 



fc=0 



where K[q N } is just the probability generating function of the binomial variable Z . This is well known 
to be 

X N , f-, XN\\ nN , , ^Xn^ ,1n 



E 



Qn 



Averaging over the initial conditions Aq with |Aq | = fcjv, we obtain 

\V N 
1 n 



E 



In 



S 



1 + (q N - 1) 



N 



N 



o(l). 



In the same way, using E[|Y" t \]/N — > 0, we get 



E 



1*0 AV t 

In 



By duality, 



E 



l + ( gjv -l))J 



E 



AT I \ 1*0 



1 + 



l + CflAr-l)^ 



0(1). 



o 1 



Letting A — > oo gives the desired result. 



□ 



For the binomial approximation, it was necessary to assume that Y /N — > 0. We now give a result 



for the case that both A t and Y t scale with N. 



Proposition 1.6. Let (X^), (Y t N ) be Markov processes with state space En that are qN -dual for some 
qN such that limjv^oo N(qN — 1) = —A G (— oo, 0]. Choose exchangeable initial conditions X^ , Yj^ G E 



N 



and suppose that both processes stay exchangeable at t > 0. Assume that the process 1 ' 1 converges 

w& 
to 



weakly to some process Y t , that N converges weakly to X t . Then (X t ) and (Y t ) are dual with respect 



H(x,y) = e 



-\xy 



Proof. We have 



E 



In 



= E 
= E 

N 



In 



{ 1 + ^r^) 



-> E 



since by exchangeability, -^t^ W^cfW — > XtYo- 



□ 
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Remark 1.7. Note that for these results we only assume duality of the processes, and not necessarily 
strong pathwise duality in the sense of However, if we start from a graphical representation, this 
result provides a pathwise construction of the limiting moment duality; see Section [3] for a related 
discussion. 

In the remainder of the paper, we discuss in some detail the case of annihilating duals and possible 
dual mechanisms. We restate Theorem 11.51 in this particular case, and discuss several examples where 
this result can be applied to rederive certain dualities, mostly known in the literature. The examples 
illuminate in particular the connection between annihilating duals and population models with balancing 
selection, as studied for example in BEM07 . In the last section we consider different values of q. The 
last example answers an open question of |AH07j concerning a self-duality derived in |AS05j . 

2 Annihilating duality 

2.1 Annihilating dual mechanisms 

In this section, we discuss the special case of a (/-duality with q = — 1, which is an annihilating duality. 
This duality relation can be written as 

P*(|Xf A Y N \ is odd) = ¥ y (\X^ A Y t N \ is odd) (5) 

for all i,j/£ En. In order to apply our rescaling result, we identify some basic mechanisms which lead 
to annihilating dualities. It is interesting to compare them to some of the coalescing mechanisms. In the 
following table, we give a list of the mechanisms that we are interested in, and afterwards discuss their 
duality relations. 





/(0,0) 


/(o,i) 


/(1,0) 






f H 


(0,0) 


(0,0) 


(1,1) 


(1,1) 


resampling 


f c 


(0,0) 


(0,1) 


(0,1) 


(0,1) 


walk-coalescence 


f A 


(0,0) 


(0,1) 


(0,1) 


(0,0) 


walk-annihilation 


f D 


(0,0) 


(0,0) 


(0,1) 


(0,1) 


death-walk 


jBC 


(0,0) 


(0,1) 


(1,1) 


(1,1) 


branching-coalescence 


f BA 


(0,0) 


(0,1) 


(1,1) 


(1,0) 


branching-annihilation 



The names given to the mechanisms are chosen to suggest an interpretation. In the resampling mecha- 
nism, the first position gives its type (0 or 1) to the second one. In the death- walk mechanism, a particle 
in the second position dies, after which a particle in the first position walks to the second position. 
Walk mechanisms suggest that a particle in the first position jumps to the second position, and either 
coalesces or annihilates if there is a particle present. In branching mechanisms, a particle in the first 
position produces a new particle in the second position, which either coalesces or annihilates with a 
particle already present. 

Remark 2.1 (Coalescing duals). In |AH07j . the coalescing dual mechanisms were classified. Concerning 
the dualities given in the above table, the following coalescing dualities were established: (i) f R and / 
are coalescing duals, (ii) f D is a coalescing self-dual, and (hi) f BC is a coalescing self-dual. 

Lemma 2.2. (a) Two basic mechanisms f,g are annihilating dual mechanisms if and only if 

\x A (g(y^))^ | is odd O \ f(x) A y\ is odd. 

(b ) With the notation of the above table, we have the following: 

(i) f R and f A are annihilating duals 
(ii) f D is an annihilating self-dual 
(Hi) f BA is an annihilating self-dual. 
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Proof, (a) is obvious. We verify (b) using the table of basic mechanisms, (i) We have that f R (x) A y 
is odd if and only if x = (1,0) or x = (1, 1), and y = (0, 1) or (1,0). In both cases, (f A (y r )Y = (1,0), 
and (1,0) A x is odd if and only if x £ {(1,0), (1,1)}. By (a) this proves the duality of f R and f A . 
For (ii) note that f D (x) A y is odd if and only if x £ {(1, 0), (1, 1)} and y £ {(0, 1), (1, 1)}. But then 
(/ D (2/ t )) t = (1,0), and the claim follows, (Hi) For f BA (x) A y to be odd we need x = (0,1) and 
y £ {(0,1), (1,1)}, ori= (1,0) and y £ {(0, 1), (1, 0)}, or x = (1,1) and y £ {(1, 0), (1, 1)}. In all cases 
(/ BA (y'l'))t a a; is odd, and there are no other possibilities. □ 

Remark 2.3. The list of duals is not complete. For a full classification of coalescing duals see |AH07] . Note 
that f R and f D have both a coalescing and an annihilating dual process. The death-walk-mechanism 
f D is g-self-dual for any gei: From the table of dual mechanisms one can check that \x A (f D (y^))^\ = 
\f D (x)Ay\ for all x,y£{Q,l} 2 . 

Remark 2.4. It is easy to see that g-dual mechanisms always satisfy /(0, 0) = (0, 0). However, unlike the 
case of coalescing duality, a mechanism need not be monotone in order to have an annihilating dual, as 
can be seen from the self-duality of the branching-annihilating mechanism. 

We can now restate our Theorem 11.51 in the special case of annihilating duals. This special case is 
motivated by the observation, made in |BEM07] , that a particular model of populations with balancing 
selection, after a transformation of the form x 1 — 2x, is dual to a double-branching annihilating 
process. Our result shows why this transformation occurs in annihilating processes. A non-spatial 
version of this model will be discussed as an example in the following section. 

Corollary 2.5. Let (X^), (Y t N ) be Markov processes with state space E N such thatF x (\X? AY N \ is odd) = 
¥ y (\X^ AY t N \ is odd) holds for all x,y £ En- Let fc^r, un £ N, and choose exchangeable initial conditions 
xn,Vn £ En such that \xn\ = fcj\r, \vn\ = n N- Suppose that both processes stay exchangeable at t > 0, 
and assume that un/N — > and E[|Yj^|/A?] — > as N — > oo. Then 

(i-^f\ 

provided that the limits exist. 

2| X N | 

As before, assuming that a limiting process (pt) of 1 and n t of \Y t ^\ exists, these processes 

satisfy the moment duality 

Proof. Corollary [53] is a consequence of Theorem II .51 by setting q = — 1. It can also be understood from 
the fact that the probability that a binomial random variable with parameters n,p takes an odd value is 
given by ^ (1 — (1 — 2p) n ) . Then we have, by binomial approximation, 

F(\X t N A Y N \ is odd | \X?\ = x N , \Y N \ = n N ) = l(l-(l-^ + o( i), 

as in the proof of Theorem 11.51 again duality, averaging over the exchangeable initial conditions, and 
taking limits, gives the result. □ 

2.2 Examples 

In this section we derive some (mostly well-known) dualities by rescaling dualities of interacting particle 
systems. We will assume that the following mechanisms occur in the process (X^) : f R occurs with rate 
^ for each ordered pair (i,j), i,j £ {1, N}, f c with rate f A with rate 2 N L , f D with rate ~cl, 

f BA with rate and f BC with rate §k Moreover, set b N := b a N + b° N . 

Consider the process l-X^I taking values in {0, ...,N}. Note that if \X^\ = k, then the number of ordered 
pairs of certain types is easily computed: The number of (0, l)-pairs (or equivalently of (l,0)-pairs) is 



lim E 

N-yoo 



2\X< 



ivixi^r 



N 



lim E 



G 



equal to k(N — k), the number of (1, 1)— pairs is equal to k(k — 1). Hence, the process \X^\,t > 0, makes 
the following transitions: 

k^k + 1 at rate rN ^ N k(N - k), (6) 
k^k-1 at rate ^±^-k(N - k) + ° N + d » + k(k 1), (7) 

k^k-2 at rate jrH k ~ !)• ( 8 ) 

In the next sections, we will consider processes of this type and their duals for various values and scalings 
of the rates. 

2.2.1 Branching annihilating process 

Let qn = djv = b c N = 0, and assume ^ — ► a > and 6^ = 6^ — > /3 > 0, as AT — > oo. The different 

scaling of the mechanism is interpreted in the sense that in the limit, the resampling affects pairs of 

\x N \ 

particles, while branching happens at a fixed rate per single particle. The rescaled discrete process £ 
has, according to © and ([7J), discrete generator 

&/ (*)_^_*,(,(^) + ,(i-)_ v (* 



A" \ \ JV / \N J J N v ' V \ iV / J 

Assume now -4 — > x as A" — > oo and / twice diffcrcntiable. Then, noting limAr^oo A" (/ (-^r-) — / (77)) 
/'(*) and limjv^^ iV 2 (/ (*±1) + / (*^) - 2/ (£)) = .f'(x), we see that G N f(k/N) converges to 

G/(a;) = /3z(l - 2x)f(x) + ax(l - 
which is the generator of the one-dimensional diffusion given by the SDE 



dX t = 0X t (l - 2X t )dt + y/2aX t (l-X t )dB t . 

This is a Wright-Fisher diffusion with local drift (3x(l — 2x). The drift has the effect of pushing X t 
towards the values and 1/2 and may be interpreted as a selection promoting heterozygosity - this 
interpretation will become more evident in the next example. Note- that it is not difficult to incorporate 
death as well: If — > 5 > 0, the resulting diffusion reads 



dX t = 0X t (l - X t )dt - SX t dt + y/aX t (l - XtjdB t 

a r"»t~\ art a o 4~ ro+o 

N 



Consider now the dual process. According to Lemma [2.21 (Y t N ) where f A happens at rate r ^-,f BA at 



^ is an annihilating dual of (X^ 1 ). The generator of \Y t N \ is 

G N f(k) J-^k(N - k) (f(k + 1) - f(k)) + h fk{k 1) (f(k 1) - f(k)) 
+ r JL k (k-l)(f(k-2)-f(k)). 

As A~ — > 00, when f(n) — > fast enough as n —> 00, this converges to 

Gf(k) := 0k (f(k + 1) - f(k)) + ak(k - 1) (f(k - 2) - f(k)) , 

which is the generator of a branching annihilating process on No- Including death, we get 

Gf(k) := 0k (f(k + 1) - /(*)) + ak(k - 1) (f(k - 2) - f(k)) + 5k (f(k - 1) - f(k)) . 

Weak convergence of the rescaled processes is obtained from the convergence of the generators by classical 
theory, see for example [EKj . By Corollary [23] we obtain for the limiting processes (Xt), (Y t ) the duality 

E x [(1 - 2X t ) Y °] = E y [(1 - 2X ) Y <] . 



Remark 2.6. Write pt := 1 — 2A" t . Ito's formula yields dpt = 0(pt —pt)dt — y a(l — pf)dB t , from which - 
at least heuristically - it is easy to read off the moment duality of the process (pt)t>o and the branching 
annihilating process. Our method establishes this duality in a straightforward manner, and also shows 
why the transformation p t = 1 — 2X t has to be applied. 
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2.2.2 Double-branching annihilating process and populations with balancing selection 



One of our motivations was to understand the transformation x i— > 1 — 2x applied in [BEM07 in order 
to obtain the duality between the competing species model and double-branching annihilating random 
walk, which is parity preserving. This is substantially different from our last example, as a branching 
event produces two new particles and not one, which is not taken care of by our basic mechanisms. 
However, it is easily implemented if we allow for multiple arrows in the graphical representation, or, 
equivalently, for basic mechanisms / : {0, l} 3 — > {0, l} 3 . 

Assume that for each ordered pair the J' 4 — mechanism happens at rate and construct an 
additional mechanism / in the following way: For each ordered triple k), i,j, k = 1, N, draw, at 
rate an arrow from i to j and from i to k. Then, if an arrow is encountered, a transition f BA occurs 
for the two pairs and (i, k). This means that at such a double transition, the state of the triple 

(xi,Xj,Xk) is changed according to the following table: 



X 


(000) 


(001) 


(010) 


(100) 


(101) 


(110) 


(011) 


(111) 




(000) 


(001) 


(010) 


(111) 


(110) 


(101) 


(011) 


(100) 



The dual mechanism / is given by inverting the arrows and applying the dual mechanism f BA to each 
of the two arrows, that is, to the pairs (xj,Xi) and (xk,Xi) with the additional rule that two l's at the 
same place annihilate each other, that is, given by the table 



X 


(000) 


(001) 


(010) 


(100) 


(101) 


(110) 


(011) 


(111) 


m 


(000) 


(101) 


(110) 


(100) 


(001) 


(010) 


(011) 


(111) 



It is clear that these two mechanisms are annihilating dual mechanisms, either by direct verification, or 
by noting that the double-branching transition is the result of two f BA — transitions happening one right 
after the other, cf. Figure 2 for a graphical representation. 




Figure 2 

Let now (Y t N ) be the process constructed from the graphical representation, where f A happens at 
rate and / at rate Then \Y t N \ has the transitions 



k + 2 at rate 
k — 2 at rate 



k(N -k)(N 



1), 



N 2 

^k(k-l)(k-2) + ^k{k-l), 



since k(N — k)(N — fc — 1) is the number of (100)— triples if there are k l's, etc. Assume ^ — >• a and 
6at— ^/3asiV— s-oo. Then the generator of \Y t N \ converges to 

Gf(k) = pk(f(k + 2) - f(k)) + ak(k - l)(/(fc - 2) - f(k)), 

which is the generator of a double-branching annihilating process. For the dual process (X^), with 
mechanisms f R and /, we obtain the transitions 

k^k + 1 at rate ^2k(N - k)(N - k - 1) + ^-k(N - k), 

fc->fc-l at rate ^2k(k - 1)(N - k) + ^-k(N - k), 
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which yield for TV — > oo, if — > x, 

Gf(x) =20x(l - x)(l - 2x)f'(x) + ax(l - x)f"(x). 

G is exactly the generator of the non-spatial version of the competing species model of |BEM07| , given 
by the SDE 

dX t = 20X t (l - X t )(l - 2X t )dt + y/2aX t (l - X t )dB t . 

For a motivation of this model as well as results on the long-term behaviour of its spatial version, see 
[BEM07] . 

3 Other values of q and self-duality of the resampling-selection 
process 

Alkemper and Hutzenthaler |AH07j ask whether the self-duality derived in [AS05] for the so-called 
resampling-selection process 

dX t = sX t (l - X t )dt - mX t dt + v / 2rX t (l - X t )dB t (9) 

could be constructed pathwise using the approach of dual basic mechanisms. This question can be split 
into several distinct questions. First, is it possible to construct (X t ) as the scaling limit of self-dual 
processes (Xf) for interacting particle systems in such a way that the self-duality of (X t ) is inherited 
from the self-duality of (X^ ) (compare |Sw06| )? Second, is it possible to explain the self-duality of the 
discrete process (X^) using a pathwise construction? Third, can we choose to construct the discrete 
processes with q-dual basic mechanisms, thus obtaining interacting particle systems that are strongly 
pathwise dual? 

As we shall see, the answer to the first two questions is yes: There is a pathwise construction, 
using the basic mechanisms from Section [21 yielding g-dual processes (X^) and (Y t N ) that rescale to 
resampling-selection processes. In general, however, our basic mechanisms will not be g-dual in the sense 
of Dcfinition ll.il the third question remains, therefore, open. 

The g-duality of our interacting particle systems will not follow from g-duality of basic mechanisms 
but, instead, from a criterion derived in |SL95j . applied in a context similar to ours in |Sw06j . Let us 
first explain how this criterion applies to processes constructed from our basic mechanisms. Sudbury 
and Lloyd consider Markov processes on {0, 1} A , for some graph A, whose generator is of the form 

Gf(x) = Y,q(i,j)(^x(i)x(j)(f(x -Si-Sj)- /Or)) + - x(j))(f(x + 6,)- f(xj) 

+ cx{i)x{3){f{x - <y - f(x)) +dx(i)(l - x(j))(f(x - Si) - /(*)) 
+ ex(i)(l - x(j))(f(x ~5 l + 6 3 )~ /(a))), x 6 {0, 1} A 

with non-negative parameters a, . . . ,e, and q(i,j) defined as follows. When i and j are neighbors in A 
(meaning that they are connected by an edge in the graph), then q(i,j) = 1/TVj, with TVj the number 
of neighbors of i; when i and j are not neighbors, q(i,j) — 0. Thus when A is the complete graph on 
{0,1,..., A}, q{i,j) = 1/TV for all i ^ j. The letters a, &, c, d and e refer to annihilation, branching, 
coalescence, death and exclusion. 

We are interested in Markov processes with state space {0, 1}{°<---' N 1 constructed from the basic 
mechanisms of Section [5] and rates chosen as follows: For every pair the mechanisms f A , f BA , 

f BC , f C , f D and f R happen at rates a/N, b a /N, b c /N, c/N, d/N and r/TV. The infinitesimal generator 
of this process is of the Sudbury-Lloyd form ([3]) with 

a = 2a, b=b a + b c + r, c = b a + c + d, d = d + r, e = a + c+d (10) 

and q(i,j) = 1/TV for all i ^ j. We shall refer to this process as the process obtained from the basic 
mechanisms via the rate parameters a, b a , b c , d and r. Note that not every Sudbury-Lloyd process can 
be constructed with our basic mechanisms: for example, if 2e < a, any solution of Eq. (|10j) has negative 
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rate parameters c < or d < 0. Furthermore, the construction is not unique - there are in general 
several choices of rate parameters a, . . . , r leading to the same Sudbury-Lloyd process. 

Sudbury and Lloyd give several conditions for g-duality of their models. A concise formula is [SOOl 
Equation (9)] , which in our notation reads 



a' = a + 2q7, & =^ + 7i c' = c — (1 + 0)7, d = d + 7, e' = e — 7 (11) 

where 7 = (a + c — d + bq)/(l — q). Eq. (fTTj) is easily translated into a criterion for processes obtained 
from our basic mechanisms. 

Proposition 3.1. Let (X t ), (Y t ) be the Sudbury-Lloyd processes obtained from our basic mechanisms 
with respective rate parameters a, b a , b c , c, d, r and a' , b a ' , b c ' , c', d' , r' . Then 

1. (X t ) and (Y t ) are dual with parameter g € R \ {1} if and only if 

a = a + <ry, b a ' - r' = b a - r, b c ' = b c + 7, c' - r' = c - r - (2 + 9)7, d' + r' = d + r + 7, 
where 7 = (2a + (1 + q)b a + qb c + c - (1 - g)r)/(l - q). 

2. (Xt) is self- dual with parameter q if and only if q = (r — 2a — b a — c)/ (b a + b c + r) . 

Proof. Eq. ([TT]) gives a necessary and sufficient condition for g-duality of two processes on the same graph 
A with infinitesimal generators G and G of the form ([3]) , see |S00] . The rates a, . . . , d for the generator 
G are obtained from the rate parameters a, . . . , r of the basic mechanisms via Eq. (|10[) ; similarly for the 
primed variables. Plugging the expression of a in terms of a, etc., into Eq. (| 1 1 [) proves the first claim. 
For the second claim, we note that (X t ) is self-dual if and only if 7 — 0, which determines q as a function 
of a, . . . , r. □ 

Now we turn to the resampling-selection process. Fix r, s, m > 0, set ajy = c^ = b°^ = 0, and choose 
rates b c N -> s > 0, aV -> m > 0, and r w /A ->• r > as N ->• 00. Fix < > and let (A^ v ) < r <t, (i / J - Ar )o<r<t 
be processes constructed from our basic mechanisms with the same rates b c Nl aV and r^. As explained 
before Definition 11.11 we can construct these processes on a common probability space, using the same 
driving Poisson processes; (X^) uses them in the forward time direction, (Y^) uses them backwards. 
The discrete rescaled processes (|A^ v |/A^)o< r <t and (|y r Ar |/A)o< T <t both have formal generator 

0»j(±)='-&«N-V(f( k -^) + f( k -^)-lf(t 



N ) N \ \ N J \ N J \N 



(I. 



For N — > 00 and fc/A^ -> x, this converges to 

Gf(x) = rx(l — x)f"(x) + sx(l — x)f (x) — mxf'(x), 

which is the generator of the diffusion ([9]), and one can show that the rescaled processes converge to 
two resampling-selection processes (X T ) and (Y T ). Now, by Proposition 13.11 (X^) and (Y^) are dual 
with respect to qN = r N /{b c N + rjv), and by our assumptions on the rates, we have N(q^ — 1) — > —s/r. 
Proposition 11.61 therefore yields 



E, 



e -(s/r)X t y 



Ey 



e -(s/r)xY t 



which is the self-duality of the resampling-selection process proven in |AS05j . Thus we have provided a 
pathwise construction of this self-duality, and we have shown that the self-duality can be obtained by 
rescaling dual interacting particle systems. 

We should note that the latter fact was shown by Swart |Sw06j . His argument, however, starts from 
independent discrete processes (X^ ) and (Y t N ), see the proof of Proposition 6 in |Sw06| . In contrast, our 
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pathwise construction provides a non-trivial coupling of (stopped) underlying discrete processes, which 
might be of interest in some contexts. 

To conclude, we mention that the duality relation can be understood from the following computation. 
For simplicity, we drop the iV-dependence in the notation. Suppose that at a given time r, there is an 
arrow from i to j, but we do not know of which type it is. The duality parameter q = r/(b + r) is nothing 
else but the (conditional) probability that this arrow is of the resampling type. We may think of the 
arrow as a randomized mechanism with transition matrix 

P(x,x') = ql x , =fR{x) + (1 - q)l x , =f Bc {x) , 

and the self-duality is explained by the relation 

Vx, y e {o,i} 2 : Yl P(^'W x ' Avl = E p (y,y'W xAy ' 1 - 

x'e{o,i} 2 y'e{o,i} 2 
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